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In this article, we study the finite temperature Casimir effect for scalar field with Robin boundary 
conditions on two parallel plates in a background spacetime that has a compact internal manifold 
with arbitrary geometry. The finite temperature Casimir force acting on the plates is derived using 
the piston approach, with which we show that the divergences of the Casimir forces from the region 
between the plates and the region outside the plates cancel each other. The sign and asymptotic 
behaviors of the Casimir force at different limits such as small and large plate separation, low and 
high temperature, are studied in detail. The range of Robin coefficients where the Casimir force 
are always attractive or always repulsive are determined. There are also some range of the Robin 
coefficients where the Casimir force is shown to be attractive at small plate separation and repulsive 
at large plate separation, or vise versa. 
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^ ; I- INTRODUCTION 

^: n 

Although Casimir effect has attracted a lots of attention due to its important roles in various fields of physics [1[ , 
there have been relatively few works on Casimir effect for quantum fields subject to Robin boundary conditions. Robin 
boundary conditions arise naturally in the study of Casimir effect for electromagnetic fields in perfectly conducting 
£S) , spheres [2j. However, in parallel plate configuration, perfectly conducting boundary conditions only lead to Dirichlet 
' boundary conditions and Neumann boundary conditions. Nevertheless, it was shown that in some geometries, Robin 
boundary conditions can be considered as extensions of perfectly conducting boundary conditions, where the Robin 
coefficients are related to the skin-depth parameter describing the finite penetration of the field into the boundary 
@; El @ ■ Robin type conditions are needed for conformally invariant field theories in the presence of boundaries, 
• ■ since Robin boundary conditions can be made conformally invariant, but ordinary Neumann boundary conditions 
cannot. The Casimir effect for fields with Robin boundary conditions was studied from different perspectives and 
under different context in [ESS B B EH EE EH EH f° r parallel plate geometry, in [l4| for a D-dimensional ball, 
in [lj| for spherical shell geometry, and in [161 ] for the configuration of two spheres. On the other hand, quantum 
fields with nonlocal boundary conditions can also be reduced to one with Robin boundary conditions [17]. The 
importance of the Robin boundary conditions to spacetime models and quantum gravity has been emphasized in 
[la , EH. For Randall-Sundrum spacetime model [2(| Hi]], Robin boundary conditions are natural for scalar fields 
and fermion fields, where the Robin coefficients are related to the curvature scale and the boundary mass terms of 
the fields [12, HH 0, H3L_ Casimir effect in braneworld models with Robin boundary conditions on the branes was 
discussed extensively inlillMllilllllSIMIMISIMIMIMIlSlll. 

Despite the importance of the Robin boundary conditions in field theory, the research on Casimir effect with Robin 
boundary conditions has been limited to a small community, and very few of the works has taken into account the 
influence of finite temperature. The goal of this paper is to consider the thermal Casimir effect for the parallel plate 
configuration. Since spacetimes with extra dimensions have become prevalent in physics, we are going to discuss 
along this line. Our setup is the same as in [l3[, i.e., we consider a pair of codimension one parallel plates in the 
background spacetime M dl+1 x Af n , where M dl+1 is the (di + l)-dimensional Minskowski spacetime, and J\f n is a 
compact internal space. In fact, Casimir effect on parallel plates in a background spacetime of this form has been 
investigated in [H, [3^, H3, [4^, [H, [44[ , when the boundary conditions are either Dirichlet boundary conditions 
or Neumann boundary conditions in the case of scalar fields, and perfect conductor boundary conditions in the case 
of electromagnetic fields. In the cases considered in [H, [H, [Hj], the dimension of the macroscopic space is 

d\ = 3 and the internal space is assumed to be the simplest one-dimensional compact space S 1 or the more general 
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FIG. 1: A movable piston inside a closed cylinder divides the cylinder into two chambers. 

n-dimensional toroidal manifold T n . In [43l l44j. the internal space is generalized to arbitrary compact manifold and 
the parallel plates are allowed to be finite with arbitrary cross section. The temperature corrections to the Casimir 
effect in the scenario of [ill . [44| were considered in [4^, 0, [47| . It has been shown that for massless scalar field with 
Dirichlet boundary conditions on both plates or Neumann boundary conditions on both plates, the Casimir force 
is always attractive, at any temperature and for any spacetime geometry. For massless scalar field with Dirichlet 
boundary condition on one plate and Neumann boundary condition on another plate, the Casimir force is always 
repulsive. The present work can be considered as a generalization of (l3j to take into account the finite temperature 
correction, or a generalization of [47| where the boundary conditions are the more general Robin boundary conditions 
instead of the Dirichlet boundary conditions or Neumann boundary conditions. We are mainly interested in studying 
the dependence of the Casimir force on the interplay between the spacetime geometry, temperature and boundary 
conditions. 

As in pil . l46l l47j , we first consider the Casimir force acting on a piston which moves freely inside a closed cylinder 
with arbitrary cross section (see FIG.[1|. In Section UH the Casimir energy in the two regions divided by piston is 
computed using cut-off method. In Section IIII1 the Casimir force acting on the piston is derived. The limit where 
one end of the cylinder is moved to infinity is equivalent to two parallel plates embedded orthogonally inside an 
infinitely long cylinder. Letting the cross section of the plates become infinitely large give the usual infinite parallel 
plate configuration. This piston scenario introduced by Cavalcanti [48| has the advantage that the contribution to the 
Casimir force due to the vacuum fluctuations of the quantum fields in the region between the plates and the region 
outside the plates are both taken into account, and the divergences of the Casimir forces from the two regions will 
cancel. In Section |IV[ we study in detail the sign of the Casimir force at small and large plate separations. 

In this article, we use the units where h = c = fee = 1- 



II. THE CASIMIR ENERGY 

As in [H, H3|, we consider a massive scalar field <p(x) in a background spacetime of the form M dl+l x Af n , where 
M dl+1 is the (d± + l)-dimensional Minkowski spacetime and J\f n is an n-dimensional compact internal manifold. Let 

ds 2 =g^dx^dx u = r] a pdx a dx 13 - G ab dy a dy b , 

/1, v =0, l,...,d; a,j3 = 0, 1, . . . , d\) a, b = 1, . . . , n, 

be the spacetime metric, where rj a p = diag (1, — 1, . . . , — 1), y a = x dl+a for a = 1, . . . , n and G a bdy a dy b is a Riemannian 
metric on Af n , and d = d\ + n. The equation of motion for a massive scalar field f(x) is 

[ J—d^yJteWdv + mA <p(x) = 0. 



We will first consider the Casimir force acting on a movable piston inside a closed cylinder of arbitrary cross section 
as shown in FIG. [TJ Let the cylinder be the region [0, Li] x SI x TV" in R dl+n , where fl is the cross section of 
the piston, which we assume to be a simply connected region in R dl_1 . The position of the piston is denoted by 
x 1 = a. If the right chamber of the cylinder is infinitely long, i.e., L\ — > 00, we obtain the configuration equivalent 
to two codimension one parallel plates located at x 1 — and x 1 — a embedded in an infinitely long cylinder. For 



3 



the boundary conditions, we assume that the field ip(x) satisfies Dirichlet boundary conditions on the curved surface 
[0, Li] x dfl x JV n of the cylinder. The interesting part is the boundary conditions on the three plates perpendicular 
to the x 1 -direction, which we take to be Robin boundary conditions: 
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= 0, (1) 



where oti,at2, 0:3 are nonnegative dimensionless constants, and ft, ft, ft are constants with dimension of length. For 
i = 1, 2, 3, a, and ft cannot be zero simultaneously. ft = corresponds to the Dirichlet boundary conditions, ai = 
corresponds to the Neumann boundary conditions, and the general Robin boundary conditions only depend on the 
ratio f3i/ot.i. We would like to remark that our ft/a* is equal to the — ft in [131 ]. 

To find the Casimir force acting on the piston, one does not have to calculate the Casimir energy of the region 
outside the cylinder, since it is independent of the position of the piston 48]. Therefore, we only have to compute the 
Casimir energies inside the left chamber and the right chamber. It is easy to see that the Casimir energy inside the 
right chamber can be obtained from the Casimir energy inside the left chamber by replacing a with L\ — a, ol\ with 
013 and ft with ft. 

Using separation of variables, one can show that the eigenmodes of the field f{x) confined in the left chamber 
satisfying the boundary conditions specified above are given by 

<Pkj,l( x ) = ^ {Asinzkx 1 + Bcosz k x 1 ) faj (x 2 , . . .,x dl ) 4>mM- 3 G N, I G N =NU{0}. (2) 

For j — 1,2,...,, (j>Qj(x 2 , . . . , x dl ) is an eigenfunction with eigenvalue j > for the Laplace operator with 
Dirichlet boundary conditions on dft. For / = 0,1,2,..., </>j\f,i(y) is an eigenfunction with eigenvalue ujJj- , for the 
Laplace operator on J\f n . By convention, when I = 0, <?W.o(2/) is a constant function with eigenvalue luJj- = 0. The 
boundary conditions JTJ) imposed on x = and x 1 = a imply that 

-f3 lZk A + ai B = 0, 
(a 2 sin(azfc) + ftz fc cos(az k )) A + (a 2 cos(az k ) - ftz fc sin(azfc)) B = 0. 

This system gives a nontrivial solution for (A, B) if and only if z k satisfies the following condition: 

F{z) = (a\OL2 — ft ft .z 2 ) sin az + (a^ft + aift )z cos az = 0. (4) 

The zeros of F(z) are either real or purely imaginary. As is shown in 0], these zeros are all simple. In the following, 
we assume that F(z) does not have imaginary zeros. This is the case if [3]: 

{ai=0,/3 2 >0} or {a 2 = 0, ft > 0} or { ai a 2 ^ 0, ft > 0, ft > 0} or [ ai a 2 ^ 0, — + — + a < 0, ftft < 1 . 

I Oil a-2 J 

(5) 

Notice that z = is always a solution of ([4]). However, it is easy to check that z k = gives to a nontrivial function 
AsvazkX 1 + B cos z k x x satisfying ([3|) if and only if a\ — 012 — 0, i.e., in the case of Neumann boundary conditions 
on both plates. We let < Z\ < z 2 < z 3 < . . . be all the positive solutions of (01 . Then the set of eigenfrequencies 
{^k.j.i} for the field tp(x) is given by 



uk, 3 ,i = J4 + w nj + w lf,i + m2 = V z l + m li > k,jeN, leN , (6) 



where 



m li = ufi^ + uj 2 m j + m 2 . 



In the case ot\ = 012 = 0, we have to let k starts from instead of 1, and zq = 0. 
The finite temperature Casimir energy inside the left chamber is given by 

OO OO OO S -. N 

^as(A) = E E E \^e- X —< +T\o g (l - e -^' T ) , (7) 
k=i j=i 1=0 ^ ' 

where A is a cut-off parameter. As was shown in (4(| up to the term of order A , 

p l m >F^ T (d + 1 - i) ^ log[A/z]-^(l)-log2 + l j- T ( i no .r i( 2^ m \ 

^Ca S (A) = > ~ C cyM A H — = C cy i, d+ i - — (Ccyl,T(, U J + lo g[M JCcyl,r(UJJ , (8J 
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where [i is a normalization constant with dimension length 1 , £ cy i t{s) is the finite temperature zeta function defined 
by 



oo oo oo 



Ci,rW = EEE E (^+m|, + [2 7 r P T] 2 )- s , 

fc— 1 j — 1 /— p— — oo 

and c cy u, < £ < d + 1, are heat kernel coefficients defined by 

oo oo oo c?+l 

EEE e-*w-^.o - E + (*) as * - 0+ - ( 9 ) 

k=l j=l 1=0 i=0 

It can be shown that (see Appendix [C]) these coefficients are linear functions of a. Moreover, the coefficient of a in 
Ccyi,z is independent of the Robin coefficients (Xi, fit, i — 1,2. 

Using a more general form of the generalized Abel-Plana formula [H, [5(| (see Appendix |A|) . we compute Ccyi,r(0) 
and Ccyir(O) m Appendix 151 The results substituted into ([8|) give the Casimir energy in the left chamber: 

E^ S (X) =S£(A)+aSi(A) 

' (A + [27rpT] 2 - ai) (ft ^/m 2 z + [2^T] 2 - a a ) 

1 — r r ; r-i 



CK. 



vEE E ^ i- ^ : = = (rkri = = = ^ v<^ 
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j=l 1=0 p= — oo 



(pi^m^ + ^-KpTf + a-^ (/? 2 ^m 2 ; + [2^T] 2 + a 2 ) 



(10) 



where Sq (A) and Si (A) are terms that are independent of a, and Si (A) is also independent of the Robin coefficients. 
In fact, Si (A) can be interpreted as the vacuum energy per unit length in the x 1 direction that would present in the 
region between the plates if the plates are absent. 

Our assumption that F(z) does not have imaginary zeros will make each term under the logarithm in (|10p nonzero. 
However, if f3\a\ < or /3 2 a 2 < 0, we may get terms that are infinite. To avoid this kind of complications, from now 
on we only consider the case where (3\ > and /3 2 > 0. 

III. THE CASIMIR FORCE 

For the piston system shown in FIG.[U the Casimir force acting on the piston is given by 

F^ oa (a, Li; a ,0) = (^ as (A) + 25&-(A)) . 

Upon differentiation with respect to a, the terms Sq (A) and S^(A) that are independent of a will be killed. On the 
other hand, for the term proportional to a in E^ as (X) + E^ as (X), we find that it is equal to aSi + (L\ — a)Si = LjEi, 
which is independent of a. Therefore, this term will also be killed after differentiation with respect to a. Consequently, 
all the terms that would diverge when A — > + will not contribute to the Casimir force. Therefore, we can set A = 
after taking derivative with respect to a to obtain 

Fg£ m (a > Li;a,/3) = - A (££ M (A) + E%J\)) = F^ s (a; a x , (h, a 2) (3 2 ) ~ F<? as (a, L r , as, fo, a 2 , fa), (11) 
where 



00 00 00 \/ m i I + [2vrpT] 2 
Fi(a;a 1 ,P 1 ,a 2 ,P 2 ) = -T V V V - t , V , ' r ^ = ^ = . (12) 

3 P °° (ft v 7 " 1 " i + l 27r P T l 2 - Q 1 ) (ft V m j, I + l 27r P T l 2 " Q 2 ) 

and 

- F ci s (« I L i; Q! 3,/?3,a 2 ,/3 2 ) = ^cas( £ i - a;«3,/?3,a 2 ,^2)- 

Notice that even though different Robin conditions are imposed on the walls x 1 = and x 1 = L\, the divergent terms 
of the Casimir energy will not contribute to the Casimir force acting on the piston. Therefore the Casimir force acting 
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on the piston is independent of the regularization procedure, as in the piston system for scalar fields with Dirichlet 
or Neumann boundary conditions. A crucial reason why this is so is that the term Hi(A) is independent of the Robin 
coefficients. 

It is easy to see that FQ as (a,L 1,0.3, (33,0.2, fa) vanishes as L\ — > 00. In other words, Fc as (a; a\, (3\, 02, P2) can be 
interpreted as the Casimir force acting between two parallel plates embedded orthogonally inside an infinitely long 
cylinder. In the following, we will drop the superscript L when we discuss the Casimir force between parallel plates. 
Setting (3\ — fi% — or o\ — 012 = or (3-y = 012 = in (fT2")l . one obtains respectively the Casimir force acting on a 
pair of parallel plates with Dirichlet-Dirichlet , Neumann-Neumann or Dirichlet-Neumann boundary conditions. The 
results coincide with the results derived in [47j . 

Notice that the Casimir force between parallel plates Fc as (a; a.\, (3\, a.2, (32) (H2) comes from the derivative with 
respect to the third term 



A ^as = - £ £ £ log 1 - r, =5 = = v e V 



-[27rpT] 2 



in (HU]), which was called the interaction term [Tj]. Usually this is also regarded as the renormalized Casimir energy 
between the plates. It has the property that it vanishes when the plate separation a goes to infinity. 

The expression for the Casimir force between parallel plates (fT2|) shows that the Casimir force is always attractive 
(negative) if the same boundary conditions (i.e., (3i/ai — ^2/0-2) are imposed on both the plates. This is a special case 
of the theorem [5l|, [HJ which states that the Casimir force between two bodies with the same property is attractive. 
From (fT2|) . we also find that the Casimir force decays exponentially as the separation of the plates a becomes large. 
As a function of the temperature T, (| 1 2 j) shows that the high temperature leading term of the Casimir force is linear 
in T, given by the sum of the terms with p = 0. On the other hand, let r = (Vo^fi)) 1 /^ 1 ^ 1 ) and R = (Vol(A/"")) 1/n be 
the size of the cross section fi and the size of the internal space N n respectively. Recall that m 2 ; =w^.+ cjjy l + m 2 . 
Since wn,j oc 1/r and ujfj^i oc l/R, (|12p shows that the Casimir force -Fbas(a; ot\, Pi, 012, P2) between a pair of parallel 
plates goes to zero when the size r of the cross section goes to zero or the mass m goes to infinity. When the size R of 
the internal space goes to zero, all the terms with I ^ go to zero, and the limit of the Casimir force is the Casimir 
force in the (d\ + l)-dimcnsional Minkowski spacetime given by the sum of the terms with 1 = 0. 

For the low temperature behavior, we use the Abel-Plana summation formula (|A3j) with 



/(z) = -2T££ 



[2ttTz] 



\ ^ (pl^rnl l + [2KT^+ ai )( l 3 2 ^ m l l + [2irTz\*+ a2 ) 2a^m? , + \2ttTz]* _ ' 
(fh V 'ml l + [2^T Z ^~a 1 )(f3 2v / m l l + [27TTz]^-a 2 ) 



The term 



OC: 







f(x)dx 

on the right-hand side of (|A3j) gives the zero temperature Casimir force acting on the parallel plates: 



1 x 



\j x2 + m li 



oo 00 „oo 

\. l.t 

dx 



7T ~{^Jo {Pi V^ 2 + m f7+°i)(feV^+™|£+Q2) c 2ajx 2 +rnj , 1 

x 2 dx 



oo oo /.oo 



7T ^ — — ' / {l3ix+ ai )(f3 2 x+a 2 ) 2ax 1 !~2 2~ 



The term 



gives 



OC: 



I o2ny _ I y 



;EE/ m ^f^-lEEE^m 



F i to m i-' F i to p=l 
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This term is independent of a and the Robin coefficients aj,ft, i — f , 2. Finally, it is easy to check that the poles of 
f(z) are exactly at 



with residues 



Z = ± ^W Z * +m £ | ' = l,2,...,i = 1,2,. ..,Z = 0,1,2,..., 



7 2 

Res '■ k 



2tt 



/li 7 -2 ( ' , aift , «2fe A 



Therefore, 

^ Res z=i j / /(z) - Res z= _,; y /(z) ,°° 

7TI 



p 2ny _ 1 — 2-^i 2-^i 



" 27TV - 1 ,,r,,, I7 2 + rr> 2 (a 4- "t& + .. ^ Z' tM+™T 

J = 1 fc=1 V ^ V + Phl+ci) ex P ( T 

From these, we find that the temperature correction to the Casimir force is 

rj-, OO OO OO 

7 \ \ \ ^ m 3,l ts fP m 3l 



A T F Cas (a; a u ft, a 2 , ft) = - - £E E — *i ( 



p V T 

j=l z=o p=l ^ 



OO OO OO 



z 2 1 ( 14 ) 

j=i z=o fc=i yjz k +m jtl [a+ pizi+ai + p , zi+ai ) exp f v ^ j -1 

In the case of di = 3 with Dirichlct-Dirichlet boundary conditions, ft = ft = and = irk /a, we find that this 
formula agrees with the formula (3.7) we derived in [461 ]. From (|14[) . it is obvious that the temperature correction 
goes to zero exponentially fast when the temperature T approaches zero. When a <C 1/T, since zu oc 1/a, the leading 
term of the temperature correction is the first term on the right hand side of (fT4|) that is independent of a. 

The next thing we would like to investigate is the limit of the Casimir force when the cross section of the plates 
is infinitely large, i.e. Vol(fi) — > oo. In this case, we have to consider the Casimir force density acting on the plates, 
which is defined as 

•Fcas(a;ai,ft,a2,P2] = am . 

r— >oo VOl(S2j 

This limit is given by replacing U)a,j by u> and turning the summation over j G N into an integral: 

vskg^ ^v^rm / (15) 

After a change of variables, we find that the finite temperature Casimir force density acting on a pair of infinite 
parallel plates is given by 



di -3 

OO OO 



T ^ ^ f°° (x 2 -m 2 ~[2TrpT} 2 ) 2 x 2 dx 

J- Ca ,(a;a 1 ,ft,a 2 ,ft) = - ^ ^ nft .^) ch . ~ ' (16) 



where 



The corresponding interaction term of the Casimir energy density is 

rp OO OO ,.00 



A£ Cas (a; ai ,ft,a 2 ,ft)=- di _ x - V V / (x 2 - m 2 - [2^pT] 2 ) 2 z 



(ft a; + ai)(ft?x + a 2 ) 



(17) 
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One can check that in the case of Dirichlet-Dirichlet, Neumann-Neumann or Dirichlet-Neumann boundary conditions, 
(HHJ) agrees with the results derived in [471 ]. 

Applying the prescription (|15p to the zero temperature Casimir force (|13|) . we find that the zero temperature 
Casimir force density acting on a pair of infinite parallel plates is 



-T=0/ 

dx 



•^Cas (a;ai,/?i,a 2 ,/? 2 ) 

, . I y " ii r 

zdlO 



t / w dl-2 / X 

V 2 / (=0 V 1 (pix— a 1 )(f3 2 x — a 2 ) 



"•' - 1 v^ 2 ~uj 2 -mf Qg\ 



1 r * 2 (* 2 ~m 2 )^ 

This agrees with the result of [l3j . The finite temperature Casimir force density is the sum of the zero temperature 
Casimir force density (fT8| and the thermal correction Ar^Cas given by 



Ar^Casiasai,/?!,^,^) = d +1 > > / l -K 1 \ v 1- 

oo oo » m 2 i 

2 Z k 1 



-i ^ ^ poo 

A yy / w * 



+ w 2 + to 2 f a + «A4hr + i ) p™ f 
V k 1 \ Pi z l+ a i f 3 2 z k+ a 2 J ex P 1 T 



-duj 



i=0 p=l 

di -2 



2"Tr-7r' li 2— \ /' / - \ T ) 



l oo oo oo o / AT — I 2\ 2 / , 



(19) 

In the massless case, the sum of the terms with I = in the first term on the right hand side of (|19|) has to be replaced 
by 

We observe that for massive case, the temperature correction term to the Casimir force density acting on a pair of 
infinite parallel plates is exponentially suppressed. However, for the massless case, the leading term of the thermal 
correction is of order T dl+1 , and this leading term is independent of the boundary conditions imposed on the plates. 

The low temperature expansion of the interaction term of the Casimir energy density can be obtained by directly 
integrating (fT8j) and (fl~9f . It is given by 



A£ C as(a;ai,/?i,a2,/3 2 ) =— — 37 — 7— Y\ \ x(x 2 - m 2 )^ log \l 



[fixx — a\){P2X — a 2 ) _ 



ii±l 00 00 , N i±i 
2 2 7T 2 J= =1 V P J 



e 2a2; 



^-,^1 00 00 00 / y — 9 — ~~ 

T 2 / ^ + 



The leading behavior of the finite temperature Casimir force acting on a pair of parallel plates when R <C a -C r 
(i.e., the plate separation is much larger than the size of the extra dimensions, but much smaller than the size of the 
cross section) can be obtained from the corresponding leading behavior of the Casimir force density acting on infinite 
parallel plates. We need to consider the case of high temperature and the case of low temperature separately. In the 
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high temperature regime, aT ^> 1. The leading term of the Casimir force when JJ«a«r and am 1 is obtained 
from the I = p = term in (116|) : 

d 1 -3 

TVol(fi) f°° (x 2 -m 2 ) 2 x 2 dx 



F Cas (a) ~ - 



2dl _ 2 il T l T fd^lX J m (ftx+ ttl )(ftx+ n2 ) 2aI _ j 

V 2 / (^l£C-ai)(^2X— «2j 



di -3 



TVol(fi) f°° (a; 2 - (am) 2 ) 2 x 2 dx 



2 <i 1 -2 7r ^i r (di_i) adl J am (f x+ai )(f x+a2 ) _ ■ 

If (3\ > and /32 > (both non-Dirichlct conditions), then in the limit a <C (i\ and a <C 02, we find that the leading 
term is 



TVol(fi) f°° x^dx _ (d x - l)r (f ) Cn(di) TVol(O) 



2di-2 7r ^j— p ^Ol_Lj a dl e 1 2 dl 7T~ a 

This leading term is the same as for the case of (3\ = 02 = (Dirichlet-Dirichlet boundary conditions). If f3± = and 
/?2 > (one Dirichlet and one non-Dirichlet) , then in the limit a <C /?2, the leading term is 

rvoi(n) /- 00 ar"'- 1 ^ _ (d 1 - i)r (f) c*(di) r ol _ dn rvoi(n) 

^ w ~ 2rfl - 27r ^ r( ^i )adI y e^+i 2 , 17r ^ ^ J *- (21) 

Notice that if both plates have non-Dirichlet boundary conditions, the leading behavior is the same as both plates 
having Dirichlet boundary conditions. In this case, the Casimir force is attractive when R <C a <C r. On the other 
hand, if one plate assumes Dirichlet boundary condition and the other assumes non-Dirichlet boundary condition, 
then the leading behavior is the same as the Dirichlct-Neumann case, i.e., the Casimir force is repulsive in the limit 
R < a < r. 

In the low temperature regime, i.e. aT -C 1, the leading behavior of the Casimir force when R <C a -C r can be 
obtained analogously from f) 1 8[) . We find that when both plates assume non-Dirichlet boundary conditions or when 
both plates assume Dirichlet boundary conditions, the leading term is 



F cas(«) 2dl+ljT ±±± ^TTT- (22) 



When one plate assumes Dirichlet boundary condition and the other plate assumes non-Dirichlet boundary condition, 
the leading term is 

^. ^ffP ,,-^! ,23, 

These zero temperature asymptotic behaviors have been observed in [l3l |. 

The asymptotics (I2TJ1) . (|2"TT) . (f2"2")l and are derived under the assumption that iiCa, i.e. the size of the extra 
dimensions are much smaller than the plate separation. In the case R ~ a, one has to take into account the correction 
terms from the I ^ terms in (p~6|) and (fl8| . In the other extreme where a <C i?, the extra dimensions play the same 
role as the cross section in the macroscopic spacetime. Therefore, when a <C R, the asymptotics of the Casimir force 
are obtained from (|20|) . (f2Tjl . ([22| and ([23| by replacing d\ with c? = di + n. We see that the sign of the Casimir force 
is not changed when we pass from R -C a to a -C i?, as long aso<r. 

At first sight, it might be quite surprising to find that when the plate separation is small, the leading behavior of 
the Casimir force when non-Dirichlet boundary conditions are imposed on both plates is the same as when Dirichlet 
boundary conditions are imposed on both plates; but the leading behavior of the Casimir force when Dirichlet boundary 
condition is imposed on one plate and non-Dirichlet boundary condition is imposed on the other plate is the same 
as when Dirichlet boundary condition is imposed on one plate and Neumann boundary condition is imposed on the 
other. In fact, this can be explained as follows. The behavior of the Casimir force with respect to a is governed by 
the solutions Zk of the function F(z) ([J|, i.e., zt satisfies 

e 2taz k = (ai ~ iPiZk){a2 - ifozk) 
(ai +z/3iz fc )(a 2 +i/3 2 Zfc)' 
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When k — > oo, Zk — > oo. Therefore if 0\ > and 02 > 0, 

e 2laZk ~ 1 as k -> oo. 

Consequently, we find that when fe is large enough, 

7r/c 
a 

is close to the corresponding for the case where 0\ = 02 = 0. On the other hand, if /?i = and 02 > 0, then when 
k is large enough, 

e 2iaz k J as fc _> oo. 

This implies that 

7f(fc-l/2) 

a 

is close to the corresponding Zk for the case of Dirichlet-Neumann boundary conditions. 

IV. ANALYSIS OF THE SIGN OF THE CASIMIR FORCE AND ITS APPLICATIONS 

In this section, we analyze in more detail the sign of the Casimir force. The cases of Dirichlet-Dirichlet (0\ =02 = 0), 
Neumann- Neumann (a\ = ai = 0) or Dirichlet-Neumann (0i = 0, oti = 0) boundary conditions have been studied 
in [47| . where it was proved that for Dirichlet-Dirichlet or Neumann-Neumann case, the Casimir force is always 
attractive. For Dirichlet-Neumann case, the Casimir force is always repulsive. In the following, we will not consider 
these cases. 

First we consider the case 0\ — 0, a-i ^ and 02 > where Dirichlet boundary condition is imposed on one of 
the plates, and generic Robin boundary conditions on the other. In this case, the finite temperature Casimir force 
between the plates is 



oo oo oo l iyy,2 



, , , -r [2npT] 2 

%(a;D;« 2 ,ft)=r^^ £ , (24) 



1 = 1 1=0 p=—oo v ; J , e V 3,1 1 r J _L y 

fe^ I + [2 I rpT]2- ct2 C 



and the zero temperature Casimir force is 



FEr s °(a;D;a 2 ,02)=lj2J2 

~ A—i r_n 



^^4e 2 ^ + 1 . /x*-m 2 ' (25) 



3,1 



If 

< g < mm ym^ + pirpT] 2 } = ^ 2 n i +m^ (26) 

then each term in the sum of (pM]) and (|^5|) is positive. We find that the Casimir force between the plates is always 
repulsive and the Casimir force is a monotonically decreasing function of a. Moreover, the presence of the extra 
dimensions enhances the Casimir force. The condition PB")) means that the ratio 0.2/ '02 cannot be too large. In other 
words, this is a close-to-Neumann condition. Therefore, it is reasonable that the Casimir force is repulsive as in the 
Dirichlet-Neumann case. In general, we have shown in Section [IIII that for any 02 > 0, the Casimir force is repulsive 
when a <r. When aT ^> 1, the leading term of the Casimir force is determined by the p = I = 0, j = 1 term in (|24|) . 
which is repulsive if ([26| is satisfied, and attractive if |26|) is not satisfied. At zero temperature, (|25|) also shows that 
if (p?6"|) is not satisfied, then when am ^> 1 oro/r> 1, the Casimir force will eventually become attractive. Therefore 
we see that if (I^B)) is not satisfied, the Casimir force will change from repulsive to attractive at any temperature. 
In the massless case, the right hand side of (|26|) goes to zero when the size of the cross section r goes to infinity. 
Therefore, for a pair of infinite parallel plates with Dirichlet boundary condition on one plate and non-Dirichlet and 
non-Neumann boundary condition on the other plate, the Casimir force always change from repulsive to attractive 
when a increases from to 00. 
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FIG. 2: The Casimir energy Ai?cas(a) and Casimir force -Fcas(fl) due to a massless (m = 0) scalar field when the macroscopic 
space is three dimension, i.e., di = 3, in the presence of the internal manifold S 1 (T 1 ) with radius R. Here the cross section of 
the plates is a square [0, L2] x [0, L3] with L2 = L3 = lm. In this figure, Pi/a 1 = 0/im, $2/0-2 = 0.3/xm and T = 0. 




FIG. 3: Same as FIG. H but with T = 1 TeV. 



An example of the first case is shown in FIG. and FIG. [3J where the graphs for the Casimir energy and Casimir 
force acting on a pair of parallel plates embedded in an infinitely long rectangular cylinder with cross section lm x lm 
due to massless scalar field is shown. The internal manifold is a circle with radius R. In these figures, we show the 
variation of the Casimir force with respect to the plate separation a at T = and T = 1 TeV (1.16 x 10 16 K), in the 
case that the internal manifold has radius R = 0.2/xm, R = 0.5/im and R = 1/Ltm and in the case without internal 
manifold (R = 0). One of the plates assumes Dirichlet boundary condition, i.e., fix = 0, and the other assumes Robin 
boundary condition with (3 2 /a2 = 0.3^im. Notice that a 2 / 'fa = 3.33 x 10 6 m _1 > wn ; i = nm^ 1 . These graphs show 
that the Casimir force is repulsive at small a and becomes attractive for a > a c . The critical point a = a c is a 
minimum point of the energy. Therefore, it is a stable equilibrium point. The position of the minimum point can be 
affected by the size of the extra dimension and temperature. The temperature dependence of the Casimir force when 
a = 0.5/^m and a = 1/xm are shown in FIG. [4] They verify the linear dependence of the Casimir force when aT > 0.5. 
Notice that with plate separation a = 1/im, the energy between the plates is of order 10 -3 TeV at T = and of order 
10 10 TeV at T = 1 TeV. The Casimir energy increases by a factor of 10 13 from T = to T = 1 TeV. Therefore the 
temperature correction is important in the high temperature regime. Experiments on Casimir effect have not been 
able to reach this high energy regime. However, with the advent of technology, one can expect that future experiments 
will be able to explore the high temperature regime which might bring forward some new applications of Casimir 
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FIG. 4: The dependence of the Casimir force on temperature when a = 0.5/im and a = 1/im. The other parameters are the 

same as in FIG. [2] 
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FIG. 5: Same as FIG. [3] but now the internal manifold is T 1 ,T 2 , T 3 , T 4 , T n is a direct product of n circles each with radius 

R — 0.2/im. The other parameters are the same as in FIG. [2] 



effect in technology. Return to FIG. IH we also notice that when the size of the extra dimension change from 0.2/im 
to 1/im, the Casimir force change from attractive to repulsive. As we have discussed in the previous section, this 
cannot happen for Dirichlet-Dirichlet, Neumann-Neumann or Dirichlet- Neumann boundary conditions. In FIG.[5l we 
show the dependence of the Casimir force on plate separation when the internal manifold is T 1 — S 1 , T 2 ,T 3 and T 4 
respectively, where T n is the product of n-circles each with radius R — 0.2/im. This figure shows that contrary to the 
Dirichlet-Dirichlet, Neumann-Neumann or Dirichlet-Neumann cases, the increase in the number of extra dimensions 
can reduce the magnitude of the Casimir force. 

The second case where a\ = 0, a.2 ^ and 02 > 0, i.e., Neumann boundary condition is imposed on one of the 
plates, and generic Robin boundary condition on the other is exactly the opposite of the first case. More precisely, in 
this case the finite temperature Casimir force and the zero temperature Casimir force are given respectively by 



\2npT} 2 

-tyy y v j " — . 

j = l 1=0 p= — oo v 7 t 1 ; = == e V 3< l 1 e ' — 1 



(27) 
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and 



F^°(a;N;a 2 ,fa) 



1 oo oo 



dx 



(P2x+a 2 ) 2ax 



1 



(28) 



Therefore, we conclude immediately that when a 2 /fa is small enough to satisfy (f2"6")l , the Casimir force is always 
attractive as in the Neumann-Neumann case. However, if a 2 / 'fa does not satisfy (|26p . the Casimir force is attractive 
at small plate separation, but will eventually turn to repulsive when aT or a/r or am is large enough. In the case 
(|2"6")) is satisfied, we can say more. As in the Dirichlet-Dirichlet or Neumann-Neumann case, (|2T))) implies that the 
magnitude of the Casimir force is a monotonically decreasing function of the plate separation. Moreover, since each 
term in the sum of (|2"T|) and (|2"5|) is positive, and the sum of the I = terms in (|2"T|) and (|2"5|) corresponds to the 
Casimir forces without extra dimensions, we find that when ((26)) is satisfied, the presence of extra dimensions enhances 
the magnitude of the Casimir force. 
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FIG. 6: Same as FIG.0 but with /3i/ai = 0.4m, f3 2 /a 2 = 0.3/xm and T = ITeV. 




Finally, we consider the generic case where a\ > 0,fa > and a 2 > 0, fa > 0. Without loss of generality, assume 
that fa/a 2 < fa/ai. As is already observed in Section HTT1 when the Robin conditions on the two plates are the same, 
i.e., faja.\ = fa/a2, then the Casimir force is always attractive. If a 2 /fa satisfies (|2"o| . so does ai/fa. The Casimir 
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force is then always attractive and its magnitude is enhanced by the presence of extra dimensions. For any fii/ai > 
and fa /ct2 > 0, we have shown in Section lllll that the Casimir force is always attractive when the plate separation 
is small enough. In the other extreme where a/r 1, if the temperature T is not zero, the dominating term of the 
Casimir force is given by the term with p = I = and j = 1 in (jT2J) . We conclude that when a is large enough, then 
the Casimir force is repulsive if 



Pi 



< 



+ 77T < — 



a 2 



and is attractive otherwise. If the temperature T is zero, the same conclusion can be derived from (Qji 

FIG. [5] and FIG. [7] show the dependence of the Casimir force on plate separation when T = 1 TeV. In FIG. the 
internal manifold is S 1 with different radius R. In FIG. [71 the internal manifold is T 1 ,T 2 ,T 3 ,T 4 respectively. The 
Robin coefficients in these graphs are fii/ai = 0.4m and 02/02 — 0.3/im. Notice that Qi//?i = 2.5m~ 1 < ixvaT l < 
3.33 x 10 6 m _1 = a%j 0ivaT x . The graphs show that the Casimir force is attractive at small a and becomes repulsive 
at large a. There is only one equilibrium point which is unstable. 

Another interesting case is shown in FIG. [H with fii/ai = 0.08/j,m and P2/&2 = 0.3/xm. In this case, wjj.i < 
ai/fli < 0,21^2- The graph shows that the Casimir force is attractive at small and at large a as dictated by our 
analysis above. However, the Casimir force can become repulsive at some intermediate values of a. This implies that 
there are two equilibrium points a± and 0,2, one unstable and one stable. If the initial separation of the two plates 
is in the range of the two equilibrium points, i.e. a\ < a < 02, then the two plates would tend to repulse each other 
until they settle at the distance a — 0,2- 
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FIG. 8: Same as FIG. [2] but with 0x1 ax = 0.08/im, p 2 /a 2 = and T = leV. 



In the figures shown in this section, the size of the internal manifold is chosen to be between 0.2/im to 1/^m so as 
to demonstrate significant difference with the Casimir force without extra dimensions. Extra dimensions of this size 
is not physically interesting. One can show that if the size of the internal manifold R is ten times smaller than the 
plate separation, there is no significant difference between the Casimir force with or without the internal manifold. 
Therefore for physically interesting size of extra dimensions of order 10 _12 nm, it will not be detectable by the present 
experiments in Casimir effect which measures Casimir force between objects that are lOnm ~ lOOOnm apart. 

Before ending this section, we would like to remark on the original setup we consider - a piston moving freely inside 
a closed cylinder. As is shown by (fTT| . in this case, the Casimir force is the difference of the Casimir force between 
the left end of the cylinder and the piston, and the Casimir force between the right end of the cylinder and the piston. 
Since the magnitude of the Casimir force is very large when the plate separation is very small, and is very small when 
the plate separation is very large, we can deduce that when the piston is close to one end, whether it is attracted to 
or pushed away from that end only depends on the Robin coefficients on the piston and on that end. However, when 
the piston is away from both ends, then which side it will move to depends on the Robin coefficients on the piston 
and on the two ends. There are some combinations of Robin coefficients that will make the piston stay in the middle 
region. 
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V. CONCLUSION 



In this article, we studied the interplay between geometry, temperature and boundary conditions on the sign and 
magnitude of the Casimir force acting on parallel plates. We first derived the finite temperature Casimir force acting 
on a piston moving freely inside a closed cylinder due to a scalar field with Robin boundary conditions. It is shown 
that even the Robin coefficients on the two ends of the cylinder are different, the Casimir force acting on the piston is 
independent of the regularization procedure, since the divergent terms of the Casimir forces from the two chambers 
of the cylinder divided by the piston are independent of the Robin coefficients and cancel each other. By moving one 
end of the cylinder to infinity, we obtain the Casimir force acting on a pair of parallel plates. In the high temperature 
regime, the leading term of the Casimir force is linear in temperature, which shows that the Casimir force has a 
classical limit. A modified Abel-Plana summation formula is used to rewrite the Casimir force which is suitable 
for the analysis of its low temperature behavior. When the parallel plates has finite size, the Casimir force decays 
exponentially when the temperature tends to zero. In case of infinite parallel plates, the temperature correction is 
of order T dl+1 , where d\ is the dimension of the macroscopic space. Interestingly, these behaviors are independent 
of the values of the Robin parameters. The sign of the Casimir force when the plate separation is small and when 
the plate separation is large is analyzed in detail. It is found that if Dirichlct condition is imposed on one plate and 
non-Dirichlet condition is imposed on the other plate, then the Casimir force is repulsive when the plate separation 
is small enough. If non-Dirichlet conditions are imposed on both plates, the Casimir force is attractive for small 
enough plate separations. We give an explanation for these behaviors by the asymptotics of the frequencies. When 
the separation between the plates becomes sufficiently large, the sign of the Casimir force depends not only on the 
boundary conditions, but also on the geometry of the transversal dimensions. We show that for a wide range of Robin 
coefficients, the Casimir force can change from attractive to repulsive or repulsive to attractive, giving rise to unstable 
equilibrium and stable equilibrium respectively. This can be applied in nanotcchnology if Robin conditions is used to 
model the skin depths of real materials. 

As mentioned in the introduction, Robin boundary conditions arise naturally in Randall-Sundrum spacetime model. 
The results in this article is not readily transferred to the Randall-Sundrum model except for massless scalar field 
that couples conformally to scalar curvature. There is a brief discussion of this in the zero temperature case in [l3l |. 
It will be interesting to consider the finite temperature Casimir effect due to a bulk massive scalar field with general 
curvature coupling, and with general Robin boundary conditions on the branes. A special case has been considered 
in [lH where the thermodynamic energy was shown to have a minimum that might give rise to brane stabilization 
mechanism. 



APPENDIX A: THE GENERALIZED ABEL-PLANA FORMULA 



Here we present a more general Abel-Plana summation formula, which is a direct generalization of those presented 
in (49l [50(1 . If fo(z), fi{z) and /2(f) are meromorphic functions, and 



lim 



lim 



{fo(x + iY)- f 1 (x + iY)}dx = > 
{fo(x-iY)-f 2 (x-iY)}dx = Q, 



then 



^2 w (z)Rcs z f (z) - ^2 w 1 (z)Res z f 1 (z) - w 2 {z)Res z f 2 (z 



f)<Rc z<c 



b<Re z<c 
Im z>0 



b<Re z<c 
Im z<0 



(Al) 



J [foiu + W)- fi( u + W)} u _ b dy+ ~^ J {fo( u ~ ~ h( u - ^)} u _ b d y ~ ^~ J {fl(. x ) _ f2{x)}dx. 

(A2) 
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Here wq{z),w\{z) and W2(z) are weight functions denned by 
w (z) ■■ 

w 1 (z) ■■ 



1, 

1/2, 


if z e 3) , « 
ifzeSDo, ' 




6 < Re z < c} , 


1, 


if z e Di, 






1/2, 


if ze d®i\{b,c}, 




= {z : 6 < Re z < c, Imz > 0} , 


1/4, 


if z = b or c, 






1, 


if 2 e D 2 , 






1/2, 


ifzGc© 2 \{&, c}, 


£ 2 = 


= {z : b < Re z < c, Imz < 0} . 


1/4, 


if z = b or c, 







w 2 (z) ■■ 

This formula can be proved in the same way as in [49l . |50| . It is a direct consequence of the residue theorem. 
To recover the original Abel-Plana summation formula, let f(z) be a mcromorphic function, and define 



f Q (z) =f(z)± log (e™ - ) = «r/(z) (l + -J—^ 

h(z) =f(z)-?-log(e™)=mf(z). 
az 



If for all x > 0, 



lim /(a; ± iY)e" 



-2ttY 



we can apply the formula (|A2[) . which gives 

> + f» - r /m* + • r + ^w^w 

p=l J ° J ° j/>0 

^ y ^Z(£)__ 27ri y .*«./(*) 

Rez>0,Im2>0 Rez>0, ImKO 



(A3) 



Here we assume that f(z) does not have poles at z = n, n = 0, 1, 2, . . .. If /(z) is analytic in the right-half plane, then 
the last three terms that contain the residues of / in the right-half plane are identically zero. In this case, we obtain 
the original Abel-Plana summation formula. 



APPENDIX B: THE ZETA FUNCTION ( cy i, T (s) 
In this section, we want to compute the zeta function 



oo oo oo oo 



k=l j=l 1=0 p=~co 

and its derivative at s = 0. By definition, Zk, k = 1, 2, . . . are the zeros of 

F(z) — (a 1 ct2 — Pi(3 2 z 2 ) sin az + {a 2 j3\ + ai(3 2 )z cos az 
on the right half-plane. We can rewrite F(z) as 

F(z) = ^{(c*! + if3 lZ )(a 2 + i(3 2 z)e iaz - {a x - if3 lZ )(a 2 - i(3 2 z) e - taz y 

Let F (z) = F{z), 

Friz) = —(at! - ip lZ )(a 2 ~ ip 2 z)e- taz , F 2 (z) = — ( ai + if3 lZ )(a 2 + i[3 2 z)e l 
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and define 



oo oo oo 



/<W = EE E (* 2 + + [27rpT] 2 ) s -l g^(z), i = 0,1,2. 



It is easy to verify that the conditions (|A1[) are satisfied. Notice that Fi(z) has zeros at z = —ictj/(3j, j — 1,2, and 
i^O 2 ) has zeros at z = ioij/fij, j = 1, 2. For i = 0, 1, 2, since i*i(z) is a holomorphic function with simple zeros, the 
poles of fi(z) coincide with the zeros of Fi(z). Applying the generalized Abel-Plana summation formula (|A2j) . we find 
that Ccyl,r(s) can be written as a sum of three terms: 



Ccyl.T(s) = CcV( S ) + Ccyl.T(s) + &I,t(*)- 



The term j T (s) is independent of a: 



(Bl) 



1 2 
Ccyl,r( s ) = - ^CnxAA,r(s) + ^ 





a. 


- 


A. 



[2tt P T] 2 



)] Wi C0S(7TS) 

i=l fc,i£N,i£N ,peZ 

1 /"OO oo OO OO i 

eee (a^m^^rx; 



m 2 - [2^7f 



2„2 



(B2) 



The first term 



1 oo oo oo 

4 EE E Ki + pfcptn 8 ) 

j = l (=0 p= — oo 



comes from the zero of Fq(z) at z = 0. If cti — «2 = 0, we have to change the sign of this term to positive. The 
second and third terms in (|B2p come from the zeros of F\{z) and F<z(z). The weights Wi are defined so that if ai = 0, 
then Wi = 1/2; if a, > 0, > 0, then = and if oti > 0, pi < 0, then Wi = 1. The last term in (|B2[) comes from 
the a-independent part of 



2^/ (/i(aO - /atojdz. 



(B3) 



The term ^ t T (s) is proportional to a, coming from the a-dependent part of (|B3 



pOO oo oo oo 

/ EE E (x 2 +m|, + [2^T] 2 )- s da ; 



Ccyl.r( s ) 



'° j=i ;=o p=-oo 

a r( s -i) 



'2^ r( s ) 



It is interesting to note that this term is independent of the Robin coefficients aj,/3j, z = 1,2. Finally the first two 
terms on the right hand side of (|A2[) give £^ , T (s): 



oo oo oo 

W = ;EE E 



c, 3 



cyl.T 



3=1 / =0 p=-oo ^V^-i + I 2 ^] 



sm 



(7T S )( 



2 2 

z — m 



A _ (/3lZ-ai)(/3 2 Z-Q 2 ) c -2a Z \ 

V (/3iz + ai)(/3 2 z + a 2 ) / 



dz. 
(B4) 



This term goes to zero as a — > oo. From (|B4|) . we find that £c y i ^(0) = 0. Therefore Ccyl,r(0) = £c y i t(^) + Ccyi r(0) is 
linear in a. Moreover, the coefficient of a is independent of the Robin coefficients. The derivative of Ccyi T (s) at s = 
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can be easily computed from (|B4|) and we find that 
d,r(0)=(C^i 1 r)'(0) + (^ r i,r)'(0) 



oo oo oo 



' f/3 u /m 2 , + [27rpTl 2 - aA (f3 2 \ /m 2 , + \2irpTP 

EE E logli \ v . v r 2 v j.' 1 ^ J v B ^v^t 



a 2 



j=l 1=0 p= — oo 



(A^m^ + [27rpT] 2 + a 1 ) (fa^/m 2 ^ + [27rpT] 2 + a 2 ) 



(B5) 



APPENDIX C: THE HEAT KERNEL COEFFICIENTS c cyM 

In this section, we show that the heat kernel coefficients Ccyi^, 0<i<d+l,([9]) are linear functions of a. Moreover, 
the coefficients of a is independent of the Robin coefficients a>i, (3i, i = 1,2. 
From the theory of elliptic operators, we have 



EE' 



d+l 



j=l 1=0 /-ii 

Using inverse Mellin transform, we find that 



E cnxAOt 4 ^ + O (t' V as t -> H 



E< 

fc=i 



c-f-200 



— / T(s)t- s Ci(s)ds, 



where 



As in Appendix 151 we find that 



z 2 + to 2 ) S . 



fe=i 



C/(s) = Cf(*) + Cj>) + C, 3 (s), 



where 



Ci00 = -o"»" 



E' 



5 to 2 



?«i C0S(7TS) 



-i />oo 2 

+ i/ (. x 2 + m 2 )- s y 



of + (3 2 x 2 



dx. 



is independent of a, 



a r (s 



2^ r( s ) 

is proportional to a and independent of the Robin coefficients, and 



((3iz - ax)(fi2Z - a 2 ) 



1 f°° -s d 

(z 2 - m 2 ) S — log I 1 — -f -e 

n ./,„ az V (piz + a\){p2Z + a 2 ) 



2a z 



dz 



(CI) 













' - to 2 ") J f 




2 - m 2 J 


.A. 




A. 





(C2) 



is an analytic function of s on the complex plane. Therefore, all the poles of the function (i(s) come from C} ( s ) + Cj ( s ) ■ 
It is easy to see that the poles of T(s)(](s) are at s = 1/2, —1/2, —3/2, —5/2, . . . and all of them are simple poles. For 
C}(s), the poles can only comes from the term 



(x 2 + to 2 ) 



aiPi , A 1 

rdx 



_ x af + (3 2 x 2 



Oii 7rr(s) J 



t s-l e -tm 



-dxdt. 
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Now, 



' " dx = e- u e-^+^dxdu = ^ / . e '^ du = ^ / e -^d« 



1 + KX 2 J J Q 2 J y/t + UK K J^Jl 

(-iy P+i 



oo 

'7T t / \/7TK 



-C 



E 



k 12 ^j!(2j + 1) K J 
This shows that 

oo — tx 2 

e -t^ I _£ - dx 



1 + KX 2 

has an asymptotic expansion of the form X)j=o as f ^ + . It is then standard to show that the function 

r ( s ) f°° I 2 , 2\-s «iA 



■ f (x 2 + to 2 ) " —5 — ^-4 da; 



has simple poles at s = 0, —1/2, —1, —3/2, —2, —5/2, .... It follows that the function r(s)£/(s) only has simple poles 
at s — 1/2,0, —1/2, —1, —3/2, . . .. Applying residue theorem to (|C1[) . we find that 

oo oo oo oo 

£ e -t(,2+ m ») = £V-! ReSs=| _. (r( s )0( s )) = £^Res s=4 _ 4 (r( s )cK*)) +E* i ~* Res W- ( r ( s K?(*)) ■ 

fe=l i=0 i=0 i=0 

(C3) 

Using the fact that 

oo oo oo oo oo oo 

fc=lj = l/=0 fc=l j = l Z=0 

we conclude from (|C3[) that the heat kernel coefficients c cy i i are linear functions of a. Moreover, the coefficient of a 
in c cy i,i does not depend on the Robin coefficients on,(3i, i = 1,2 since C|( s ) does not. 
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